home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / dspcon.z / dspcon
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDSSSSPPPPCCCCOOOONNNN((((3333FFFF))))                                                          DDDDSSSSPPPPCCCCOOOONNNN((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DSPCON - estimate the reciprocal of the condition number (in the 1-norm)
  10.      of a real symmetric packed matrix A using the factorization A = U*D*U**T
  11.      or A = L*D*L**T computed by DSPTRF
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE DSPCON( UPLO, N, AP, IPIV, ANORM, RCOND, WORK, IWORK, INFO )
  15.  
  16.          CHARACTER      UPLO
  17.  
  18.          INTEGER        INFO, N
  19.  
  20.          DOUBLE         PRECISION ANORM, RCOND
  21.  
  22.          INTEGER        IPIV( * ), IWORK( * )
  23.  
  24.          DOUBLE         PRECISION AP( * ), WORK( * )
  25.  
  26. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  27.      DSPCON estimates the reciprocal of the condition number (in the 1-norm)
  28.      of a real symmetric packed matrix A using the factorization A = U*D*U**T
  29.      or A = L*D*L**T computed by DSPTRF.
  30.  
  31.      An estimate is obtained for norm(inv(A)), and the reciprocal of the
  32.      condition number is computed as RCOND = 1 / (ANORM * norm(inv(A))).
  33.  
  34.  
  35. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  36.      UPLO    (input) CHARACTER*1
  37.              Specifies whether the details of the factorization are stored as
  38.              an upper or lower triangular matrix.  = 'U':  Upper triangular,
  39.              form is A = U*D*U**T;
  40.              = 'L':  Lower triangular, form is A = L*D*L**T.
  41.  
  42.      N       (input) INTEGER
  43.              The order of the matrix A.  N >= 0.
  44.  
  45.      AP      (input) DOUBLE PRECISION array, dimension (N*(N+1)/2)
  46.              The block diagonal matrix D and the multipliers used to obtain
  47.              the factor U or L as computed by DSPTRF, stored as a packed
  48.              triangular matrix.
  49.  
  50.      IPIV    (input) INTEGER array, dimension (N)
  51.              Details of the interchanges and the block structure of D as
  52.              determined by DSPTRF.
  53.  
  54.      ANORM   (input) DOUBLE PRECISION
  55.              The 1-norm of the original matrix A.
  56.  
  57.      RCOND   (output) DOUBLE PRECISION
  58.              The reciprocal of the condition number of the matrix A, computed
  59.              as RCOND = 1/(ANORM * AINVNM), where AINVNM is an estimate of the
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDSSSSPPPPCCCCOOOONNNN((((3333FFFF))))                                                          DDDDSSSSPPPPCCCCOOOONNNN((((3333FFFF))))
  71.  
  72.  
  73.  
  74.              1-norm of inv(A) computed in this routine.
  75.  
  76.      WORK    (workspace) DOUBLE PRECISION array, dimension (2*N)
  77.  
  78.      IWORK    (workspace) INTEGER array, dimension (N)
  79.  
  80.      INFO    (output) INTEGER
  81.              = 0:  successful exit
  82.              < 0:  if INFO = -i, the i-th argument had an illegal value
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.